function [dof_x, dof_y] = build_dof_coordinate(dof_map, V, T, FE_Order)
% function dof_coordinate = build_dof_coordinate(dof_map, V, T, FE_Order)
%
%
n_dof = max(max(dof_map));
% n_elem = size(T,1);
dof_x = zeros(n_dof, 1);
dof_y = zeros(n_dof, 1);
[I, J, K] = indices(FE_Order);
px = (V(T(:,1),1)*I' + V(T(:,2),1)*J' + V(T(:,3),1)*K')/FE_Order;
py = (V(T(:,1),2)*I' + V(T(:,2),2)*J' + V(T(:,3),2)*K')/FE_Order;
        
dof_x(dof_map) = px';
dof_y(dof_map) = py';

end